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1- INTROD UCTION 

"5 e- ■ z> 



It is not always easy (of possible) to ^evaluate the 
'definite integral / a b f(x) dx by finding an «an tideriv.at iv^e , 
of f and computing F(b)-P(a)^ It may happen that, al-' 
though the integrand f is a simple function, there, is no 
simple function F such that V = f. .This happens, for ex- 
° ample, if f(x) = sinfx 2 ). Or there migjit.be 3, suitable 'f 
but the methods of .integration needed to determine F are 
excessively m tedious . evaluation of definite integrals 

occurs often enough in those disciplines where calculus * 
is used that this is an important problem/ and it can be 
solved by realizing that what: is desired is^ot the whole 
anti^erivative function F(x) , but just, the riumerical quan- 
tity F(b)'-F(a).- 

0 

We shall describe* two methods for obtaining numbers 
that are approximations to f(x) dx. How close *an*ap- 
proximation is needed is of course determined by the con- 
text in^which the integral .arises., but it is .important in 
any approximation method to have a way of knowing how far, 
at worst, the approximation may be from the true value o : f, 
the integral. For example, knowing that "f.Q999 is an ap'- 
prbximation .to k ' ' 



dx 



is a relatively useless piece of information compared to 
being able to say that j£ 

i ix = 1.099 + 0.0059 . ( 
h * , . « i * 

For the latter assertion allows one to say with certainty 
that 1 

; 1 .0940 < M i dx < 1 . 1058* ♦ 



while the original assertion really says' nothing about 
the value of the integral. 



Th-e methods covered in this unit are the . trapezoidal 
Rule and "Romberg's Method. Ttve Tripe^o'idai . Rule is a 
good place to begin for many reasons. The derivation' of 
-tke^ trapezordal formula- is relatively* straightforward; 9 
the formula \t^sel^f is easy to program on a computcy or 
programmable calculator; the, method is accurate enough i 
for many applications; and ttyere is often an easy way. to* 
determine the magnitude, at worst, oL the error. ' • 
Romberg fs Method provides a clever way of slightly modi- 
fying the results -£rcmr*he Trapezoid'al Rule therebv im- 
proving the accuracy by a tremendous amount.* Because 
these formulas involve considerable calculation, a 4 compu«- 
ter is usually used to perform these tasks. We discuss 
£ome of Vhe problems and pitfalls -of computer implementa- 
tion, such as ; rounclof f error, 'using computed- drawn graphs 
for illustration. . * , ] 

» 4 

2V THE TRAPEZOIDAL RULfr - THEORY 



* Recall,.that the definite integral ^ b f ( x ) dx is de- 
fined fusing 1 imi ts of .some sort - see Exercise 4) to co- 
incide (whe'n f(x> > 0 on [a,b]) with our intuitive' notion 

o 

of the area between the x-axis and the graph y. - f(x), 
and between x = a and x-.= b. The Trapezoidal Rule is 
♦fcased on the, observation that trapezoids may easily be • 
us.ed to approximate the area as followsV 1 

Fix a positive integer n and divide* the interval 
[a,b] into n equal subittterval s using points x 9 ,x,,..(^x 
witn-xj, = a and * n ± b. Note that x j =\ + ^iAx where 
Ax = (b-a)/n. Then by /'connecting the dots' 4 from (a\£(a).) 
to (x^ffx.j)) -to (x 2 '/f (x 2 ) , # and so on t<\{b,f(b)); one get 
n trapezoids, the sum of whose areas approximates 
jfa f ( x ) dx - Since* the area of the i tn trapezoid is 
J 2( f (x i . 1 ) + f(x i ))Ax, the' total trapezoidal area is: 

Ax((%f(a)+%£(x 1 )) + (%f(k 1 )+3sf(x 2 )) + .. .+ftf(x n *)W(b)))._ 



(x 1 ,f(x 1 



Ax 







1 


Ax 



J = f(x) 



b,f(b)). 



* a - =x 0' X l 



X 2 



Figure 1. 



Combining terms leads to the following definition. 
i)e fin jit ion . The n^T irqpezoidal approximation , T , fbr 
an integral l^(^) dx'is 



T n = *f 



i = l 



where x' = (b-a)/n, and the x^ 'are as^-stated above. 
As an example, consider the integral ^ 



■dx 



For T 



.AX = |j£ = Q.2 , J, ( < 

and the points x i are r, 1 . 2 , 1 . 4* . . . ; 2 . 8 , 3 . * Thus, ,\ 

v - . . ' * . « 

T =0.2 + -i— + -i- + + + 1 | = 'i 10156 

, »• 4-2 .1.2- 1.4 1.6 ' 2-.8J *- lui: > 0 

In this case we may use antideri\ra*tives to get 

f 3 i dx = In xi 3 = In 3 -/in 1 In 3 ■ 1.09861 . < 

3 
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8 



Thus the Trapezoidal Rule gives an approximation to In 3. 

Here T x J is not very accurate (the error is about T).003). 

With more work ^pne can " calculate T 10Q ,. which turns out to 

be c . - * 

. , T l0( = '1.09864 , 

with an error of about 0.00003. 

. i . - 

It is evident from the geometric nature of the Trap- 
ezoidal Rule that T n -is a perfect approximation whenever 
the integrand is linear (i:'e.„ has the £orm ax + b). For 
tnon-lipear functions the trapezoidal approximation is, in 
• general, not perfect, but we can predict how gre.at the 
error can be in terms of the size of the second deriva.- 
tive, f"(x), which* is a measure-of how Close f(x)* is to 
being linear. (If f(x) is linear'then f"(x) = 0.) 

Theorem . *I£ f"(x), exists on {a,b] anc} n is a positive 
integer, then there is a number c between a an.d b such , 
that '/ -v _ s. 



(2)' T - f b f(x) dx N '= f "ic) (b-a) 

" Ja 12 n z 



The point c" in this theorem^ is Qnspecified. - After 
all, if we knew c then the formula above would, allow us ' 
to cal-Culate the exact value of the. integral from T . ' 
There is no general way of telling where between a and b 
c might lie. However, the following corollary is ah im- 
mediate consequence and 'has a form that is easier to 
apply. 

Corollary . { Theoretical- error bound' fj?r T n l If\f"(x) 
exists on [a,b] and K is chosen* so; that |"f"(c)| 4. K when- 
ever a < c .< b, then " - - - N 
* ■* 

(3) ' IT - f b f(x) dxl 1 K (b-a) 3 

f I n J a . _ . I ' 1? n 2 

The proof of the theorem above is relafively long", 
and the interested] reader may consult [DM, p. 233 or 
p. 305] . , ~ * K 4 



Note, that this corollary gifes.only an upper bouivd 
on the error: Usually (although see Exercisers) the er- 
ror is below this bound, sometimes by a, considerable 
amoun t, ( , ^ t . 

Consider how the theoretical error bound applies to 
J^i'dx, In. this" case f"(x) is 2/x 3 , . which^i s posit i ve . 
In* fact, 2/x 3 is *a decreasing function on [1,5] (sketch 
its graph)wand so its maximum value tor this^interva^ \ 
occurs at x * 1, and^is .2/1 3 or 1\ Thus we may take, L 
T = 2 and apply the corollary to" deduce that the error^ett- 
T ! n will >be ,at most 2(2) 3 /12(n) 2 , or 4/3n 2 . So the^ror 

in T roo is- at most 4/30000'or 0.000133 We ^ tasted 

above that T i00 = 1.09864 and so we may now conclude that 

(4) * 1.09864 - .00014 £ f - dx < 1.09864 + .00014 « 

ox',equivalently, that 1.09850 < In 3 < 1.0V878., (Exer- 
cise 2 shows how a bit more information can be extracted 

in' thi s case . ) * t 

\ 

~ We may use the formula fpr an error bound £p find 

out which T n we should compote to^obtain a given accuracy 

Suppose we wished to calculate j^siPx dx with an error 

* o 

of> no more thin 0-001.. (Of course it is rputine to eval- 
uate fj^is integral exactly as 2, 'but it stUl serves to ' 
illustrate tl)e method.) /Since j f ff Cx) | = |- sin x f , and-* 
since sin x, always lies between -1 aifd we may 'conclude 
that, |f"(x)| < 1 and so take £ = 1 . Then by the theoret- 
ical error bound (3), we know- T that> the error T n is*at - 
-most l»7r 3 /12n 2 . Thus if n is ►so large that , 

\ Tr 3 /12n 2 <l 0.001 , 

then | n will be wijfchin the required tolerance. In sol- 
ving t^js latter inequality for n, we find that 

\ 4 • ' 

• A - « 12n 2 > lOOOir 3 4 * ' 

n '* > /100<0ir VI 2- ' t 

. * * n . y 50.8 . 



iff 



Thus for n = 51 ^(or greater)., T n has' the- desired- accuracy . 
In. fact, T 5 1 = 1.99936, which differs from 2^ 0.00^64. 

.An important consequence^o f the presence of the n 2 
in the denominator of the^ theoretical error bound in "(3) 
is t-fyat usual ly ,( though there are exceptions J doubling n. 
the^number of trapezoids, cuts the error by a factor of 
,4. Similarl^ multiplying n hy^in sJiould cut the error 
tw*a factor of 100. *For ; * 



/ 



3 1 "a 

x d \> ' / 



the. error in w T 10 is about 0.005. ifte error in T 100 is 

about O.OOOaS. ■ * 

• * 

In the previous two examples, it vfras not vpr # ]j diffi- « 
4 Quit to determine a value 'for K. The- determina t ion is ' ; 
some^tijnes a bit more ■ intricate , and on« may, have, ta use 
the methods of locating maxima. an4 minima learned oin 
elementary calculus (see Exercises 7 and 8). And* the r^ 
* are times -when the theoretical rfrrffr^bound cannot be used. 
For. v instance, if f(x) is complicate^ f"(x) may ^be even 
more so, and obtaining tfn/u^per. bound'k' for* *| f"f x) | oVi 
[a,b] may be extremely difficult. For example , 'consider 

f 8 x 3 . 

an integral that arises in thermodynamics. Here f"(x) js 
^.very complicated (chect for' yourself ) an^ a value tor K 
*is difficult to^obtain. In such a situation you might^ * e 
wish to compute manyvalues of -T^, halting when:t;hey seem 
to^have stopped changing in the^decimal plnces y.oircare 
about. This is discussed, further in the next section..' . 

i For some i^tegrMs, the error bound can be lowered by 
reducing the loss sustained'in passing from f"(c)^in 
dqua^ion (2)^ to K in the i-nequality (3). For if it is . ' 
^possible to divide [a,b] into two subintervals so that on 
one- o£*them f" (s substantially less than' the bodnd, K, i 



11 

\ . r 



for all of [a,b], then error bounds for the subijitervals 
may be computed, and added. As an example, consider 
*J 3 3? ( * x# ^ e seconc * derivative f"(x} s 2/x 3 is larger 
near 1 than it is near 3. If n is even, we may consider 
T n as the sum of for J 2 i dx and for J 3 i dx. 
Now, if Kj, X 2 denote bounds on 2/x 3 on [1,2] and [2*3] 
respectively, then ^ 

,Kj = 2/1 3 =* 2, K 2 < = 2/2 3 = 0 . 25 a 

and the -error in T n is no greater than ■ 

• . Ki I 3 A K 2 I 3 0.75 * : 
12 (Jin) * TTTW = "n*""' 

Thus the error in Tj 00 -is at most 0.000075 and the bound 
on the integral giveifl-n (4) can be improved to 

1.09857 < i < 1.09872 . / 
J j x 

The method of Exercise 2 improves this to 

1,09857 ^ ( i dx < 1.09864 . 
J j 

Another context in which the trapezoidal rule is 
useful, but the theoretical error bound less, so, is when 
the "integrand is described by a table of values, as qp- 
*t>ps ad -to a mathematical express ioit.* For- instance, values 
of a .fund: ion , f, might' be obtained in a laboratory ex- 
periment and a definite integral*-'©'/ f is required. See 
IRK] |pay^^Qf which are in [TF, p. 216] ) # o-r^[DM, p. 266] 
s for speed f ic '^applications , tjtie former* to the problem p£ 
measuring cardiac n output , the© latter -to a, ptob*leiri,/in 
thermodynamics . For example ^ sjtippos.e the following five 

values of^ A furttrtion f(x) were determined* 7 in an experi- 

i ' . - f i ■ j** 

ment, and an approximation to J 0 f(X) dx is required. ^ «• 





X 


oioa 


.25 


.50 


.75 1.0 






0.000 


.235^^ 


.388 


.420 .349 



Since .the Trapezoidal Rule needs values only at equally 
spaced points^ we can immediately compute 

,T„ = 0* 25 ^ 4, °^ 549 4- 0.235 + 0:'388^0.420 = 0.304 . 

Note tttat since the integrand is not 'given by an equation,* 
'no estimate K on the second derivative is available, and 
so the theoretical error bound is not applicable. < Be- 
cause we have no idea of the behavior of the function be- 
tween the tabulated points, > we haye no way of tending how 
close T u = 0.304 is to the true value of the integral. 
In practice one should obtain en<5ugjj data paints to elim- 
ma.te the possibility of bizarre cha~nges in the function 
.between the points. Note that t in this example, Tj = 
0.175 and T 2 = 0.381.'. In Section. 4 we shall ; see how t 
these values of T x , T 2 and T l| may be combined to yield 
,a closer approximation to the integral. 



Exercises- ' - * . * 

1. Compute T^ and T 8 for (a) J* dx; ^ (b) cos x 




- / ■ r 

-2\^(Jonsider the use of fr^ to approximate J 3 i dx. -Why is it that 
( for every n, the value 9f T n is larger than the true value of 
]' the integral? This fact implies that tne inequality in (4) can 
be. improved to 1.09850 < fn3< 1.09864. 

3. Give ah example (by sketching a graph if you like) of a function 
f and two integers m- and n with the following property: n > m 
but T^-is closer. to,,J\/f(x) dx than T n is.^ 

4. The va£ue of a definite integral can also 'be approximated by 

i 4 "Sums* of areas of xectajjgles. Fix n, and let x ,...,x partition 
[a,b] as oh page-3> Let 4 the" "sum of the areas af the n rectan- 
- gles be denoted by * * * 
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8» 



, . fn-1 • 
. R n = Hi-V^ 



(Figure 2 illustrates the case n =.4*.) 



f(x) 




R 4 = ixf(x Q ) + Axf( Xl ) + Axf(x 2 -) j Axf(x 3 ) 

Figure 2. t 

Show that T n = R n + (b-a)(f(b)-f(a))/2n 

Compute T for f x 2 dx. -How much less than, the error bound 
r t 1 * o 
given by (3) is the error in Tj? 

6- Which T n « should one compute when approximating J e^* 2 ) dx to be 
sure that the error is at most 0.0905? Show that a much smaller 
n will^do by computing the error bound separately on [0, iy and 
■ llh, 2], as was done in the" text for I — dx. 

■7. Use the inequality in (3) to show that 



1" 

J 0 



100 



If you have a computer available, calculate T to shpw that 



18.631 < [ 
" Jo 



0 x-^ 
e 6 dx <.18.662 . 



Consider the integral of Exercise 1(c). Use the inequality (3) to 
say by how. much, at worst, T a differs from the true value of the « 
integral (Hint: 'in (3), K need -not be the exact maximum of 
|f"(x)| on [a,b] (which in this case is hard to find):^ any 
upper bound on [f"(x)| will do.)., 



.14 



9. The temperature outdoors was taken every 3 hours during a 
24-hour period, with the following results: 



Time . 


mid-.- ' . - mid- 
night 3 6 9 noon 3.6 9 night 


Temp. , 


10.0° 9.1° 12. A 6 18.6° 25.9° 32.7° 31.5° 20.0° 18.9° 



Recalling that the average value of a function f(t) defined on 
[0,24] is'ij f(t) dt, use the trapezoidal rule to approximate 
5 the average 'temperature during the day. * 



3. THE TRAPEZOIDAL RULE - PRACTICE 



Because of the many calculations required to compute 
T n' °" e n % ormall 7 uses a computing machine - either a 
large scale digital computer, ox a programmable .calcula- 
tor. IF all that 'ia-wanted is T n for, a single value of ,n, 
it is very easy to write a program to compute it. How- 
ever, it is often desirable to compute many different 
values of T . instance if, as'we saw in Section 2, 

a theoretical error bound is not 'readily -obtainable , one 

would like to produce a sequence of values of T until 

1 * n 

the values cease changing by much. If one adopts* a- — 

straightforward approach to computing T/, , T T f 

n * 2 9 ioo 

the/i (2+3+. . .+101) = 5150 evaluations o£ f(x) will be 

made (since T n requires n + 1 evaluations' of f) , many of ; 

them more than once. m j \ 

* • k - / 

A more efficient way to compute -many valu^> of T n 

with no duplication .of function evaluations is^to Compute 

T j 9 T 2 ■ T «* ■ T s 9 T i 6 9 * " " ' To see wh y this is efficient 
note- that having just computed T % , say, the values. of 
f(x.) needed to compute T* are exactly the 5 values just 
used for T„ , plus the 4 values occurring midway between 
them. I # 



\' 1 -.-I 



H — ^ — I 



base points 
for T fl 



This observation" leads to the following formula for cal- 
culating T'* n from T n : " > ** > 



(5) 



'2n 



b-a 
<2n~ 



. n-1 
+ X f 
i-0 . 



a + (2i+l) 



where , the values of f at ,the evenly indexed x^ 1 s , 
^11 appear in . 



i.e. , 



Using this method one first decides the number of 
— evaluationsL._o£_£(x) thereu-is, ^ime.jEpr, say 513, and, in- 

stead of computing, just> T 5l2> or computing T , T 2 , T 3 , . . 
. T*3 0 , one can compute T 512 .and obtain the values of Tj , T 

T *» T^, T 32> T 6I>I Y l28 , T 2 5 6 along' the way: It is. 
• useful to have*all this output; for comparison. Moreover 

in Section 4 we, shall develop 'a method which begins with 

such a , sequence of values, of T n , and modifies 'them to ob 

tain still closer* approximations . 

* * Hexe is a program wrrtjen in BASIC wnich uses Equa- 



tion (5) to /generate a senuence of values of T , 
program was written fpr the integrand x 3 /(e x --l). 
other^f unctions, line 130 must be* modified, - 

100 REM TRAPEZOIDAL APPROXIMATION OF -INTEGRALS 
U0 REM FIRST DEFINE INTEGRAND; THIS LINE 
120 .REM MUST BE MPDIFIED FOR NEW INTEGRANDS • * ' 
130 DBF FNA(X)=Xf3/(EXP(X)-l) ■ 



The 
For 



ERLC 



IS 



ii 



140 DIM X(3Q) , 
ISO PEUJfT "WHAT ARE A AND B? !I 
160 I^PUT A,B 

170 PRINT "HOW MANY T(N)»S WOUW) YOU LIKE?". 
180 ifoPUT M 

190 REM ' »■ 

200 REM . WE FIRST USE, ONE TRAPEZOID, STORING 
210 REM , THE RESULT IN T(L) 

IS ^ 

220 D=B-A 

230 T(1)=D/2*(FNA(A)+FNA(B)) * 
240 REM 

250 REM NOW COMPUTE THE (2fN)TH .TRAPEZOIDAL 
,260 REM APPROXIMATION FROM THE 2f(N-l)ST, 
270 REM l USING FORMULA (5) ABO^E, AND 
280 REM STORING IT itf T^N) 
290 FOR N=2 TO M . ' 
300 D=?D/2 \ 
310 S-0 
320 REM, 

330 R^M _ USE S TO FORM THE SUMMATION ifl (5) 
340 F0R*J=1 TO 2f(N-l) STEP 2 
3S0 S=SfFNA(A+J*D) * 
360 NEXt J. ' 
370 T<(N|=T(N-l)/2+D*S 
380 NfcX^.N 

390 REMl ' 1 . 1 

REM- \ NOW PRINT THE RESULTS , ' 
410SQRi N=l to M | 
420 PRiH v ,t T( t, ;2f(N-l};' , )= l, ;T(N)^ • 
430 NEXT N 

440 END m ' 

Here i*s a sample output -(on a 12 significant 
computer)- for • # 



digit 



RUN 

WHAT ARE A AND B? • K 

HOW MANY T(N) »S WOULD ■ YOU LIKE? \ 

11 *• ' v ' 4 

.T(l)=2. 63826923395 > „ - 

*T(2)*4.9020?237702 . 

T(4)=5. 76289887395 , , 

T(8)*5.95440195p63 ' ' m 

T(16)=5. 99988421985 v 

T(32)*6.01109575704 *, ' 

T(64)=:6, 01388856817^ 

T(128)=6. 01458613933 - 

T(256}=6.pl476049262 , 

T(S12)-6*01480407847~ . , 

T(1024) =6. 01 481497477. 

i It appears from this output that 
dx = 6.0148 



to- 5 significant digits (in j^act, the true value is . 
6r.0148186« . .) * While it is generally true that, when the 
values; get close* to each other, they are close to the 
liiit, i.e.; the „true value 0f the integral, there are 
some examples pf-jtnno^ent looking integrands inhere the 
values -are] momentarily' cjuite close, but then change sig- 
nificantly as, they approach .their limit. (An example of. 
this is ^discussed in ExeTcise 17) ." 'Experience' and some 
cpnsef vatism* help* out; !e.g; , from the data above oae: ■ \ 
might; Merely conclude that' the t$ue value is /ft. 014. , , ; 
~I-£a S . significant *<te'{*it., i r s required, more values might 
be computed^ : oi another method usedV. However '1000 "is* 
often "a practical- upper bound on., the number of * trapezoids 
that can l>e used: ' 'the computer used above, working with 
so mjty , significant digits,' takes, a lot of time, while on 
a machine/with fewer significant digits the"- results are 



less trustworthy because of roundoff error, which we now 
* discuss. - ~ * ' ' 

« v • 

* «* 

So far everything that wVhav'e <Jone (foV instance, 
concluding from (3)| that T^' approaches rf(x) .dx as n + »)>^ 
is based on the fact that- .the real numbers have infinitely 
many, decimal "digits*, jfowe.ver, any.' computer can work* only * 
with a 'fixed finite n6mber of .digits $ This cireates. unex- 
pected difficulties, which of^en require great ingenuity 
to avoid. Consider ^Fi^ure which gives a pictorial . 
representation of' the results of using 'a computer* to ob- 
tain T n for j . sin x dx. .» ^ 
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2000' 3000 4000 5000, 

Number of Trapezoids ' 

Trapezoidal Error for sin (x) from 0 to 3.141592 - 



Figure 3. 



*A Digital* Equipment .VAX;, which uses 7 significant digits-, was used. 

Equation (1) was used' in this example and the one of Figure 4. 
* When Equation, (§)* is used, there^s less roundoff error, though the 

same general behaVipr ,is exhibited. * . ' » 

. •/ 'St ' , ' t * * 
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The absolute yalue^ of the error is plotted* logarithmical- 
ly on the vertical axis, against the .number of trapezoids. 

p . Id l 592 * r 

Since J 0 . sin x dx = 2.00000, this means that the ; 
graph labelled "TRAPEZOIDAL ERR0R\is the g*aph of }'T -2|\ 
The theoretical error bound for this ^ integral w.a,s worked 
out in Section 2, and is' (3 . 141592) 3 /YZn z . For small 
values ofji,.say n ^ 5O0 r the observed q\ .computed error* 
agrees quite well with the theoretical prediction - it is 
less than the^error bound, and the two graphs, are roughly 
parallel.' As is the case* with the theoretical e^ror 
bound graph, the actual error* is quarteTed when n is dou- 
bled. As the number of trapezoids, reaches and passes 
about 700, however, something quite strange, happens-. The 
actual error curve becomes less smooth and is soon e ggeat- " 
er than the theoretical error bound, in clear contradic- 
tion to*-(*3)! II one decided to compute T . , the theo- 

. < ' »-2 9 Q 0 

retical error bound would indicate an error of at most 



10 



but the error in *the computed vdXue would be 



' about 3 x 10 



- 6 



\ almost 10 times greater then expected. 
■ •** * » *>. 

, The Reason fot this deviation 'from the theoretically 
predicted behavior is that ther.computer uses "only ^sig- 
nificant digits, and the rounding necessary to .perform 
additions^in this mode can,** when repeated, Very often, r 
cajuse a substantial buildup (Of what is cdmmonl^ known* as 
roundof f'error . ... > 

For example, consider how a computer, working an f> 
significant i&gi% floating po.int arithmetic aBds',123 .4S'6 
; and .987654.* Mt';"£irst converts' to^exponential (power|qf 
ten) notation wifch'-a ctJmmqn exponent, Say .1*2'3'456.E 3 A and 
.000988 £ 3, .rounding the smaller numbe'r^it' cannot "%tore 
• 00987654 E 3. since .this requires 9 signUfican-t difgitsj*. 
Now it*adds to- 'get :\Z4aU\E 3 6'x 12t .444-; *when tjffe true ' 
sum i v s % 'i?4^ 4436*54; The eTror here (in thevam/unt of 
*.4)00346) may seem inconsequential since it does not af- : 
JEect the first*^^ignificaRtwJtigi<s of the resalty but v " - 

- .V . \9h \\ :: -cr V . - 4> ^ M . 



. wfien this occurs repeatedly tl^e error buildup Can be 

quite substantial (even though founding may cause some 

cancellation of error, since . it. 1 sometimes increases ~and f - 
• * * . * 

sometimes decreases the "fesult) . For a drastic example 

' .consider a sum of the form 1.00000 + 0.000004*+ 0.000004. 
The first. addition is done as .100000 E 1 + .000000 E 1, 
yielding a result of 1.^00000", and similarly^ the next ad- 
dition yields a final result of 1.00000. But the true 
answer is 1.00JJJ008 which, rounded to 6 digits is l.OOOOi, 
This is. exactly .what happens in the computation of T . 

* f < 5 / 1 k 1 S 9 2 ^ 

For instance, when computing T 4 0 0 for J/" -sin x dx 

^rom Equation (1) we calculate the sun/ 

399 . «" 

I sin(x . ) . 
• i=l 1 * 

The subtotal, of the first 200 terms is about 127, and in 
completing the sum we add to this subtotal 199 more num- 
' bers that are each between -1 and 1. The repeated round- 
ing causes a significant loss in accuracy. • 

' Note tRat we are using the term accuracy in a slight 
ly -new sense here, as we are not talking about the differ 

f3 1 4 I 5 9 2 * 
o * sin x dx C w hich is often 

called truncation error) , but rather the difference be-* 

-tVeen the real value of T 4(J(J and the. -computed value oi 

T^ 0 0 . It is^not> generally possible .to' determine, whether 

. Xjie roundoff error will cause the computed value of T n to 

be closer to or farther from the^actual value of the in- 

tegral. One should therefore^dUstinguish between the twjo 

types of err.or and, once it is determined that a certain 

, T wilJr have a sufficiently small truncation error, T 

. should be compyted .in a way that minimizes the roundoff 

error. We shall discuss one such technique shortly, but„, 

note that one shoulcl avoid dealing with T <at ali-M|en n 

* . a n .* ^ 

ds so large that the machine -being used- cannot caf%% 

• 6nough decimal places to calculate T n a^curatel^g^While" 

\a 13 -digit computer or programmable calcula,tox : can .com- 



pute T n where n 1000- with little loss 'of accuracy, a 
6-digit cojupater cannot. Thus if the- theoretical error 
bound says that T R is needed where n 800 (as'in Exer- 
- cise 6) and one is working -on a 6-digit machine, it is 
^ probably worthwhile- to seejc another, more efficfent 
'method of estimating the integral*, like Romberg's Method', 
•discussed in the following section, - 

Consider-again the examples cited above of a comput- 
er addition. Note that rotindoff loss is much less likely 
if the .numbers to be added' are of the same order of mag- 
nitude, ftr then the Shifting necessary to get the *expo- 
n $ n £h t0 agree is unnecessary, and the concomitant round- 
offs loss avoided.. (Rounding error is not entirely elimi- 
nated; in adding .44444 and .88888, the result, 1 .333332 
will have to be rounded to l\ 33333.) It follows that' 
When^many^umbers^w^ich vary substantially in size are to 
be added, roundoff error can be\ lessened by adding Ihem 
In ordei^from smallest -to largest . 'Th'is is because the 
smaller numbers, may then contribute to the sum and more-^ 
of- the intermediate additipns" il 1 be of two lumbers of* 
the same order of magnitude, If the larger ones were 
added first, then very few of the additions : wotfld invi^ve 
numbers of the same order of magnitude; For a simpie* ex- 
ample, note that when 1, .000004, .000004 (see above) are 
added reverse -order, " one ^first gets . 000008 , which 
is rounded to ,000001 E 1 so it can, be added tp .100000 
F}, giving- a result of 1.00001, correct to 6 significant' 
"digits^ Tn decendirig" order the computed ^sum" is 1. 00*000. 

*. This -fact is dramatically brought home" 'by considers 
ing the trapezoidal . approximation to. 'J 15 e" x dx; The - 
values o^e"» x on the* interval [0, 15] range from 1 down 
to 3 x 10r* 7 , jo substantial.rdfundoff error is expected 
for T n , the computation af which involves summing .values • 
of e-x as x ranges from U,to 15. A simple way to add- the 
numbers im reverse order, i.e., from smallest to, largest, / 



is to Consider T R for j 15 -e" x dx, where €he minus sign 
is included to account 'for the interchange of limits. 
For # the latter integral ax will be negative ,".an4 so T 
will.be computed from right to left. Theoretically T 
"for these two integrals are equal, but as demonstrated in 
Figure 4, the computed values differ by a lot when n £ ^ 
1000. Recall^that the error (vertical) scale is 4ogTT? / 
rithmic, and so* the difference is reaHy quite great. ' " 
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Trapezoidal Error for Exp(-X) ; from 0 to 15 
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'Figure 4. ' 



The mathematical- study of Voundof f error is general- 
ly quite difficult. Techniques of probability theory 
come into play, to .account for the probable cancellation 



of founding erroTs. /For instance, it can be shown that*** 
in the addition>sf n numbers, roundoff error will usually 
build up proportionately- to /n. This helps to explain 
the shape of th e ri ght hand sections of the graphs in 
Figures 3 and 4. " > 

Sometimes thfc way in which "'roundoff error creeps' 
into a*calculation is quite subtle, and anyone using a . 
computer to perform many ari thmeticaj. ^operations must be 
aware of how it can destroy certain *ypes of* calculat ions . 

• • • - 

Exercises 



IrO. In Figure. 4 why is* the theoretical error bound the same for' 
both / 0 e~ x dx and /° -e" x dx? 

H. Write a program in a language other than JSASIC, e.g., in the 
language, of a^programmable calculator, that c'ompufces^ a 
Tj ,T 2 ,T,^ ,T 8 , . . . for a given integral* 

12. Apply the program of the previous exercise to the integral of 
4 Exercise Then do it "backwards" as in the text and compare 
. the jfesults, 4 ' 



IMPROVING" THE TRAPEZOIDAL RULE BY ROMBERG 1 S METHOD . 

5 7" : 1 : 

Often rules for bounding the error in a numerical 
method are ignored in practice because they are too dif- ■ 
ficult to apply. Still, a detailed analysis of error can 
be useful^ in another way, for it may lead to a new method 
which i& better than the original one. In this section 
we will develop Romberg's Method, the basic ideas of 
which have a wide range of application in numerical 
analysis. 

Suppose a,b, and f are given, and let 
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E n - T n 1 ( b f <x>dx « f,, t cHb;a^ m 
11 " ja ^ . % 12 n\ 

As you know from Section 2, the point c depends on n. If 
We changed n to n + 1 , then most likely fe would change as 
well.- It turns out 'that E n can be represented in a dif- 
ferent, rather more 'useful way; namely, if f is suffi- 
ciently differentiable, then -* ' ' ^ 

(6) * E - V + — * ~ 

n n? n- n* 

where A, B, C, ... are constants that depend on. a,rb and" 
f , but are independent of n. "~ « 

Itl other words, ■ * * , • . 

v ' 

(7) * T n = ( b f(x)dx + * + Ll + O.. 

n ia / n 2 n\ nT \ 

This representation of the trapezoidal* er ro*\ is a conse- , \ / 
quence of the Euler-Maclaurin summation formula - see * ' 

[DR, p. 108 and p. 327], With no further/ assumptions on >i ; \ 

A, B, C,... art: seems that' the term A/n 2 will lie tie Jt-atfe-. > 
est contributor to the error. The key, idea in Romber^^**": " ^ 
Method is to eliminate this, term in thi* following wayV ^, :+ ^ % 
'We first' repl.ace n by 2n in (7) to get| . ' " - ; ^ ^V^V- 

h X„^„ I A B / C 



(8J • T^f f b ftx)dx I J± 
1 ^a 4n z 



16n^ 64n* . 
We then subtract T from 4T ? . to get ' 



o 



i„ i rum 4 I oti 

n 2n 



(9) 4T 2n - T n *- 3^f(x)dx 
When we divide both sides of Shi s~ la^'e^a^^^^5^we^ 



"*-Aq^ually; * the^xght~j>ana e^res$j&ri^^ei|^ *^?V?^ 
^"Hs^i'^^^^i^^^V"^®*^^ ^of^jig?.^aXi^ last tin & 



(10) 



* 2 



■15C/48 



+ , 



. a . n* - n° 

Equation (10) shows that * (4T 2n - T n )/3 is a better appr ox 
.amation.to the integral than either T„ or T. 9ti , in the 
sehs^ that its error has £erms involving n" fy , n" 6 , . . . , 
but'lione ..involving n~ 2 t* 

7 *If we.^let T n ' = (4T 2n -- T n )/3, it follows that the . , 
error for T£ can be written as a senies beginning with 
n"* > and we can repeat this/procedure starting with T' 
and T£ n to eliminate the n" % term from the error.. The 
• numbers 4' and 3 will be replaced by 16 and island so we 
' Ut T" = # (16T^ n r Tjp/15. ' There is no need to stop heje. 
The following diagram .shows that by -beginning, with v the. s 
*n values, T 4 i* T 2 , T % ,JT $ ,..», T 2 n-1* one .gets a jtr;La'n~gular 
array, ending with ,a row- 



7m 



T • 

1 2 



.T, 



1 1 6 



T 2 n-1 



T^ 



Tit 
*2 



(4/ 



Figure 5" 



coi^istiiig of just T^ 11 " 15 , which is usdally the/be t st ap- 
proximation in. Jhe array* ' ' \ , 

There are many reasons why this method/ is an, excel- 
lent £ne to use. Unlike many methods (Trapezoidal Rule, 
ftectang^fe Rule/ Simpson's Rule) whifch ar.e perfect ;for 
polynomials of a certain ^degree- CI » 0 v 3 •respectively), 
'Romberg's method is perfect foY*- polynomials of arbitrary 
ily. h^gh degree, i'f -one computes sufficiently- many rows. 



*More. pro&g$ely, -the first tow is perfect for 'linear func- 
tions, the second for cubics, v and in general is per- 
fect^or'. polynomials of degree c 2i + 1. ; 

Moreover, in terms- of computing effort, the amount 
of work needed tp^ apply Romberg's method once the trape- 
zoidal values' are obtained^s : really quite miitimal com- 
pared to the gain in accuracy % ; -^ndeed , if one. begins 
with^the program in Sectioif\$£ it is* a simple .matter to 
add a ffew instructions to "produce the Romberg array. No 
further evaluations of the^nte grand are needed, onjy 
some straightforward arithmetical manipulations with the 
X-(l)'s* Constructing a .suitable program for Romberg's 
k method is left as an exercise^ \ j 

Examples * 

[ 7 dx l\*3£a33 J. 16667, 1 . 11667 1 .10321 

; * - , " 1.11111 1M0000 1.09873 * 

True value: -1~.QS%26 .1.09&64 

1.09861 - ;„ 1 .,09863 . 

Vote that TpJ , Which needed T e , is closer/to the' 
true value than Tioo ('se,e, pag.e' - ) . 

* sin x dx . .0-- - 1.57080 1.89612'* 1 .97423 

V - . 2.09440 ^2^.00456 2.Q00j>7 9 

True val^e: 1.9,9857 1 1*99998 

2 • r . 2.ooooy w 

In, this case^., it take's 600 trapezoids to/ obtain a 
value for T R . as aCcurate ^s T^ 3 ? , which Required only* 
, 9 function evaluations. ** 



2.63827 4..90201 5.762T50 5.9544D 



• | „ Z,0^8Z/ 4..90Z01 5.76290 

7 ' 5.65659 6.04986 6.01324 

True value:] - 6.07608 6.0161/3 

. 6.01482 ' ? 6.0IW : '7. . 

, I-Z _Co in pari n g ^ t h the results ^in Section, 3^_we jieje^th^t;, 
here T> J \ about as accurate as T,, s . 



.This method 4 is ^particularly useful when the iiite-^ 
grand is given [by a table : of Values obtained in ah ex- . 
»peri*ient. In v the example at the end of Section 3, we 

•V - Y >;~ ■T'--- "' 27 - . . 



obtained ^ 

T, = 0.175, T =0.281, T* = 0.304 . 

1 * 2 * «♦ 

* Romberg's method may be applied to yield 

TJ = 0.316, T] * 0.312, V\ = 0.311 . 

Without redoing the experiment, that is, without obtain- 
ing more function values, we improve our or i gina 1, approx - 
imation*of th^s integral from 0.304 to 0.311. ; 

The field of approximate integration is vast, with, 
many methods suited to special types of integrands. Rom- 
berg's method is an excellent one to use when the situa- 
J tionf requires equally spaced x-values. 

Exercises / 

V 

13. Write a program in BASIC which starts with the program in 
Section 3 and produces the Romberg array. 

14. Use your program you wrote to compute the Romberg array for 

/ . T-Ldx. 

15. Compute the Romberg array for j Q cos (8 sin x-x) dx (which \ 
equals 0.73713182...). Does Romberg's method improve on the 
trapezoidal results? 

16. - If you know about Simpson* s Rule, show that the entries in the 

second row of the Romberg array (i.e., the T^) are the approx- 
' imations one gets from Simpson's Rule. j 



17. Consider the integral 

s x dx 



J. 



cosh x 



where 



1 * -x 
cosh x = j (e + e ) 



Compute Tj, T 2 , and T k . Then compute two more rows of the 
28 ■ • ' . • 23 



Romberg array. Note that TJ agrees with T" to 6 decimal 
places ^ and so one might cease the computations here, expect- 
ing the true value of the integral to be 0.479555... . But 
th'is is not so, as the true value is 0.4794282... (compute 
some more rows o¥ the Romberg array, or do the" integral by 
finding an ant iderivative) . This points out the inaccuracy 
that can arise if one stops computing when the values agree to 
a certain number of digits. (Tins example is adapted from 
[DR, p. 317] .) 

18. Compute the Romberg array for some integrals of your own 
choosing, preferably some for which you can determine the 
true value. 
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6. MODEL EXAM 



(A calculator i.s necessary for problems 1 , 2 „ 3, and 9.) 
1". With the^help of a calculator, evaluate T 5 for 
f i .i i 



e x -l 



dx. 



2. Show that if T l0 is used to approximate 



then t he > truncation error is no greater than 0.107. 

3. If you were to use T to approximate 
i 

(In x) 2 dx, 

and you wished the truncation error to be no greater than 0.0005, 
which n should you use? 

* * 

4. From the point of view of evaluating the trapezoidal rule approxi- 
mation with minimal roundoff error, which of the" following four 
forms of the definite integral . is best to use? 



a) 



c) 



i o 

MO 



dx 



dx 



b) 



d) 



x b + 1 



-2 



x* ♦ 1 



dx 



dx. 



JO 



The graph of ^ - on [-10,10] ii shown in Figure 6 



ERJC 30 




sili x dx and 

h J 



The following two integrals are equal. 

p-rr 

s in x dx. 

Jo 

lor wfjich one will the trapezoidal approximation, 1 , be 

n 

closer to the true value* of the integral? 

Why when computing a sequence of trapezoidal approximations 
to an integral, is it better to compute '] y , Tu , T u , \\ , T Jt , 
rather than Tj, T,, T 4 , 1^ , T , ...? 

T, denotes thfe first entry of the fourth row of the Romberg 
array for 
fb 



f (x) dx . 



a 




How many times must the integrand, f (x) , be evaluated in order 
to compute t( 3 ^ ° 

Which row of the Romberg array always gives the true value of 
the integral when the integrand is a polynomial of degree 5? 

SUppose all that is known about a function/f, is that f(0) = 0, 
f(l) = 1, f(2) = 3, f(3) = 5, and f(4) = 0, and that an 
approximation to 

f(x) dx ■: I 

is desired. Compute as much of the Romberg' array as is possible 
from the given set of f -values. 

Show that, when the Romberg technique is applied to I and T 

n 3n' 

rather than to T n abd T 2n , the resulting improved approximation 
has the form 

• 9 T - T 

3n n , » • ' 

8 

instead of 

4 T - T f 
2n n, , i 
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7. SOLUTIONS TO EXERCISES . 

(a) T„ = 3.13118, T 8 = 3.13899 

(b) = 0 , T 8 - 3.13899 (why'?) 

(c) T u = 1.03229, T 8 - 1.98235 

Because the graph of y - 1/x is concave up on [1,3), the 
trapezoids completely enclose more than the area under the 
graph. Thus, T n is too large. * 

















I 3 



Alternatively, f"(x) - 2/x 3 is positive on [1,3], so that the 
difference T„ - f b f (x) dx in Equation (2) is positive. 

n / a 



• 


£ X 






.51 


7i 



If f(x) as in the diagram is such that f(x) dx = 0, that is 

J o 

if the areas above and below the x-axis exactly cancel, then 
T - 0 will be perfect while T 2 is positive. 
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_ B ^ lb-a)(f(b)-f(a)) . 
" R n + 2n" 

5, r, = i( f(0) - f(1) ] = f- k ■ Since 



|f"(x)| = 2, a constant, we may take K = 2 m (3). The error 

«. 2M 3 1 
is at most — * — -r 

12 -1 6 

Since [ x^dx = ^ , the error in Tj is exactly j - ^ = g> . 
; o 

f(x) = e (x } , f(x) fc 2xe (x . } , f"(x) = (2x) (2x) fe 1 ) - 
♦ 2e (x2) = (4* a .*2W xJj . 

Since f"(x) is a positive increasing function on [0,2), 
' i f"W| <_ f M (2) - 18e* . „ 
.Using (3), we want n so that 

18e - i< ; - 23 < 0.0005, or —rr; < n 2 , or 1144.7 < n . 
12n - 0005 

Accordingly, n should be at least 1145. . 

However we may improve this estimate by noting that, for n even, 
T n is the sum of'Ti^ for 

e^ dx , 

o 

and 1, for* 

2 ? 

e^ ■'dx . 

*Since |f M (x)| < f"(1.5) = 104.365 on [0,1.5), and |f"(x)| 
<_ f"(2) = 18c 1 * = 982.^68 on tl.5,2)", the error T r is at most 

104. 365(1. 5-0) 3 + 982t 768(2-1. 5) 3 • 
12 (hn) 2 12 (hnf 

• , 158.359 * a . . ^ . 158.359 

which equals ^ . So if n is large enough to make = 

n z • n^ 

• * ' «. 

£ 0.0005, then T r will be within the iesired tolerance. 
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" - v - rT^eTgre J^su£free3-W rr?s - :r ' ' 



/ This ifer'-. ' 

eluding -ttat:;^^ ^' - ^.j^ y ' ~ r - 




and so 



f 1 ' • (x)* = a 



if 



-x 2 ^ x 2 
27 + 3" 



2X 

"3 



x 2 - 9x + 18 
(x-3)(x : 6) 
x = 3 or x 



0 

6. 



The' # maxima and minima of f"(x) must occur at x = 0^3,6, or 10± 
The' values of f" at these points are respectively, 0.667, 
-1.494, 0.667, and 0.0065, whence the maximum value of f"(-x) 
on [0,10] is 1.494, Taking K = 1.494 in (3) shows that the 



error in T 100 is at. most 

1.494 x 10 s 
% 12 x 10 2 



0.01245' 



ERLC 
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;-l'87.;,f : (x) = 



a (sin x)//I 



.-. .(cos x)e,' 

'2/2TT" 



(sin xU^2 



\ , Asm x)/yC2. r 2 - S 
- , r - W — — — ^r<v*> -sin x 



- . - Ratfer 



/2 



tb% <^archirife-for critical points of f"-(x]_, j&'jOy 
obtain; a valufe-for K-by notin^that . ^ 

.- ♦ -sin*'x] ^Tj^S|l_2iJ ^ •< ^ + J < j ^ 



and that,' oji .[0,21, 



e 



2/2-n 



2/2tj 



< 0.2283 



SoJf*(x)| < 0.2283 x 1.7072 < 0.3898 and we may let K = 0;3898 
Thus the er^or in T Q is at most 

♦ 0.3898(271) 3 

12 • 64 - 0,126 ' 



For f(t) dt, 

' n 

2£ f 10. 0^18. S 

8 \ 2 



r 

T fl 



* 9.1 + T2.4 -K..+ 20 



493.95 . 



Thus T for 



' 0 



!• is 493.95/24 = 20.6, which approximates the average temperature 
during the day. 

10. The second derivatives of the two integrands differ only in 

. f . sign. Since K must bound |f M (x)|, K is the same for both inte- 
grands . 

9 • 

11. 'Here is a program in^'the language of the T^xas Instruments 

-- SR-52- — ^..1- _ — . 

000 LBL A STO 01 HLT 
.V"006 LBL B STO 02 1 STO 06 
015 RCL 02 -^CL 01 = STO 03 
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026 X(RCL 01 E + RCL 02 E) 
038 t2= 
J^041 LBL C STO 04 HLT 

047 2 4 INV PROD 03 'PROD '6*6 
055 0 STO 05 ' ^ 

0 59 RCL 06 STO 00^ 
065 £BL D RCL 01 + RCL 03 X 
, 075 (RCL 06 + 1 • RCL 00 = 
086 E SUM 05 
090 1 -INV SUM 00 dsz D 
097 RCL 05 X RCL 03' + 
10°S RCL 04 t 2 = GTO C 
113 LBL E ... rtn 



Register usage: 


R o, 


. a 




R 02 


b 




R 0 3 


K „ b-a b-a b - a 
D -a, -g-,... 




R o. 


T T T T 

1 1 > 1 2 > 1 k * 1 8 * ' ' ' 




R 

0 s 


accumulates summation in (5) 




R nt 

0 6 


1,2,4,8,... 




R 

0 0 


index for^loop to compute 



summation 



User .Instructions - 

1* Fill-in instructions 115 on, to provide a subroutine that 
computes f(x) , the integrand. Use* register 07 if nbcessary. 
Don't use =; use (...) insteaa\ % 

Enter a, press A, enter b, press B.f *• - 

Tj is in display. Press run. , 



T 2 is in display. Press run repeatedly^to see , T Q , 
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12. Example : 
For \ 

115 STO 07 (RCL 07 - RCL 07 x' t 6)INV lnx rtn 



10 e X ' x2/6 dx, first fill in 115-132 as follows: 
o , 



Then press 0, A, 10, B, and run repeatedly to get: 

5.006363169 (Tj ) 

14.00806104' 
18.1343902 
18.^5166875 
18.61450407 
18.63898675 
18.64510821 
18.64663861 



(T l28 ) - 4 

Due to the 12-digit accuracy of programmable calculators, one 
gets the' same results for the integral backwards, i.e., 



dx. 



Discrepancies would come in with the use of thousands of 
trapezoids. On'a 6-digit machine however, the forward and 
backward results would already differ at T^ 8 . 

13. Line 140 has been changed. Lines 440 - 560 compute 
Romberg array from the trapezoidal data. 

' 100 REM TRAPEZOIDAL APPROXIMATION OF INTEGRALS 
110 REM* FIRST DEFINE INTEGRAND; THIS CINE 
120 REM MUST BE MODIFIED FOR NEW INTEGRANDS 
130 DEF FNA(X)=Xt3/(EXP(X)-H 
140 DIM T(30) 

150 PRINT "WHAT ARE A AND B?" 
160 INPUT A,B 

170 PRINT "HOW MANY ROWS OF THE ROMBERG ARRAY?" 
* 180 INPUT. M 
\ -190 REM 

s 200 REM WE FIRST USE ONE TRAPEZOID, STORING 
210 REM /'THE RESULT IN T(l) 
220 D=B^ 

230 T(1)0/2*(WIA(A)+FNA($) . 

o7 



32 



240 RBI 

250 RBI NOW COMPUTE THE (2tN)TH TRAPEZOIDAL 

260 REM APPROXIMATION FROM THE 2t(N-l)ST, 

270 REM USING FORMULA (5) ABOVE, AND 

280 REM STORING IT IN T(N) 

29(T FOR N=2 TO M 

300 D=D/2 

310 S=0 

320 REM , 
330 REM USE S TO FORM THE SlMiATIOfMN (5) 
.340 FOR J=l TO 2t(N-l) STEP 2 v ■ 
350 *S=S+FNA(A+J*D)«. 
360 NEXT* J 

370 f(N)=T(N-l)/2+D*S 
380 NEXT N 
390^ftEM* 

400 REM ^JOW PRINT THE RESULTS 

4X0 FOR N=l TO M ^ 

420 PRINT , T( ,, ;2f(N-l); ,, )= M ;T(N) ' ^ 

430 NEXT N 

^440 PRINT 

450 F=l, . . 1 

460 FOR 1=1 TO M-l 
470 F»4*F 

480 REM COMPUTE AND PRINT THE ITH ROW 
490 REM OF THE ROMBERG ARI^Y" 
' 500 ROR J=l TO M-I» 
510 T(J) = (F*^fJ+l)-T(J))/(F-l) ' 
520 PRINT T(J) ^ ' 

530 NEXT J 

S40 print 

550 NEXT I * ' * 

560 feND ' 

38 " 
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.73627360 .73713182 
.73741790 



14. 3.00000 3.10000 3.13118 3.13899 
3.13333 3.14157- 3.14154 

^.14212 3.14159 
5. 14.59 

15. 0. ()0000O0 1 3540803 .1 25062fb 
2.0721071 -.35007730 .93971077 
-.51155626 1 .0256966 . 72393171 
1.0500975 .71914179 * 

9 

.717843^2 

In this rather unusual example. Romberg's method provides 
worse approximations than the original trapezoidal approxima- 
tions. 

16. Using (5) in the form 

~ Tn , b'-a v odd 

T 2n = T + IK r - ' 



where 



5 odd v y 



l<j<2n r l 



we have 



\ - ( 4T 2n- T n^ 3 



* = (2T :> 1%^- I 0dd -T )/3 

n ,2n n J 

= I |b-a p(a)^f(b) + ^evenj + 4(b-a) ^oddj 

s W(^l( f(a) + f(b)+ ' 2even + 4 ° ad ] • / 

. which is precisely Simpson's rule with 2n subintervals . 

17. U 758664 .799332 . .559499 .499453 

.479555 .479555 '.479438 ' 

.479555 .479430 
.479428 
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1. 2.215.. 



8. SOLOTIONS TO MODBL EXAM 



2. [f"(*)| = 2/x 3 which is decreasing on [.2?iJ, so K may be taken 
to^be 2/(.2) 3 , which equals 250. Thus 

" gRROR <^L^1 = 0.10666 ...< 0.107. 

3. |f r (x)| s H I 20nx) < 2 e 613 < 2-2(lnx) < 2 a . 

4. . A * — 1 — — 



ad" 



7 



1 £ * on Taking K = ^ in the error formula yields 

n > 51.6, or n >- S2, 



4. (d) Is best suited sinqe it adds f- values in order from smallest 
*t$.lafgest ., 



nje 



5. fit * 

"a sinxdx. 



For n fixed, a smaller b --a leads to a smaller error. 



6. The first sequence reuses the f-values already competed. Moreover r 
"the first sequence can be used to start Romberg 1 s -method. 



7. 



8^ The third ^row. 

9. T x = 0 T 2 = 6 



T* = 9 



TJ = 8 



= 10 



V[ =* 10.1333 



. 10. 



T ■ 
n' 



T = 
3n 



b A 8 

f(x) dx+ 2 + 15" + ... 
n n 



t> A B 

a 4 



So 

9T - T 
sn n 



f (x) dx 



... j 



and the term of the error ^involving n 2 h^s been eliminated. 



STUDENT FORM 1 
Request for Help 



Return to: 

EDC/UMAP 

55 Chapel St. 

Newton, MA 02160 



Student: If you have trouble with fe specific part of this unit, please fill 
out this form and take it to your instructor for assistance* The information 
you give will help the author to revise the unit. 



Your Name 



Page 








0 Upper 


OR 


Section 


OR 


' ►„ 0 Middle 




Paragraph 


0 Lower 









Unit No, 



Description of Difficulty; (Please be specific) 



Model Exam * 

Problem No. 

* - 

Text, 
Problem No. 



Instructor : Please indicate your resolution of the difficulty in this box. 
Corrected errors in materials • List corrections here: 



c 



{ 



o 



Gave student better explanation, example, or procedure than in unit, 
Give brief outline of your addition here: 



( J Assisted student in acquiring general learning and problem-solving 
skills (not using examples from this unit.) 

■ . / / 
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Instructor's Signature_ 



please use reverse if necessary^ 



Return to: 

STUDENT FORM 2 EDC/UMAP 

Unit Questionnaire * f T 5 Chape i A S ^ t £fX 

x Newton, MA 02160 



Name * ' j j n ± t y 0 , Date _ 

Institution . Course No. 



Check the choice for each question t that comes closest to youx* personal opinion. 
1. How useful, was the amount of detail in the unit ? 
Not enough detail to understand *the unit 



Unit would have been clearer wieffejnore detail , 



_Appropriate amount of 'detail 



* Unit was occasionally too detailed , but this was not distracting 
Too much detail; t was often distracted 



2. How helpful were the problem answers ? 

, 

Sample solutions were too brief; I could not do the intermediate steps 

1 4 Sufficient information was given to solve the problems 

Sample solutions were too detailed; I didn't need them 



Except for fulfilling the prerequisites, how much did you use' other sources (for 
example, instructor y friends, or other books) in Qrder to understand the unit? d 

A Lot . "Somewhat A Little N ot at all ' 



A. How long was this unijt in comparison to the amount of time you generally spejxd on 
a lesson (lecture and Tiomework assignment) in a typical math or science course? 

Much Somewhat "About - Somewhat Much 
Longer Longer - ,the Same ' Shorter Shorter 



5. Were any of the following parts of the unit confusing or distracting ? (Check 
as many as apply.) * ? 

Prerequisite's + 



^Statement of skills and concepts (objectives) 
^Paragraph headings 
_Examples v * % 
Special Assistance Supplement (if- present) 
_0ther, please explain ; . 



6. Were any of the following parts of the unit particularly helpful? (Check as many 
as apply.) • ' 

^ Prerequisites 

Statement of skills and concepts (objectives) 



Examples 

* Problems 

Paragraph headings * 
_TabIe of Contents 

Sp/ecial Assistance Supplement (if present) 

Other, please explain _j 



Please describe anything in the unit that you did not particularly like. . 

• / . . / 

/ 

Please describe anything tllat^ you fdund- particularly helpful. (Please use the" back of 
this sheet if you need more space.) 



